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The probable observation of a supersolid helium phase was recently reported by Kim et al. In 
this article, we confirm their speculation. Based on a theoretical model for solid helium-4, we 
show that the emergence of such a phase is inevitable. This is the first instance of Bose-Einstein 
condensation (BEC) in a solid. We calculate the BEC transition temperature {Tc) and the Landau 
critical velocity («c) from the model, which respectively are 215 mK and 251/ims~^. They are in 
excellent agreement with the experimental results; Tc =175 mK and Vc < 300/ims~^. We also prove 
that our model possesses the necessary and sufficient condition for the emergence of supersolidity. 
We briefly comment about similar behaviour in He. 

PACS numbers: 67.80.-s, 67.90. +z, 67.40.-w 
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Helium has been one of the most exotic elements that 
has aroused an immense amount of theoretical and ex- 
perimental interest for almost a century [1|. Recently 
another possible first in the realm of Helium, the super- 
solid phase was reported |3| • This is believed to the first 
example of a Bose Einstein condensation (BEC) in solids. 
The speculations about the existence of such a phase has 
been around for more than 30 years |3j,y, [5j together 
with numerous attempts to experimentally confirm them 
i 

As is described later, helium is a quantum solid. Va- 
cancies and interstitials are essentially delocalized due to 
their ability to tunnel through potential barriers. At low 
temperatures these point defects form a weakly interact- 
ing Bose gas. Also large zero point motions result in 
exceptionally rapid exchange between helium atoms at 
nearby sites. BEC of either of these two channels may 
provide for a new phase of matter with long range crys- 
talline and supersolid order coexisting. This is the super- 
solid phase. According to Andreev, Lifshitz and Chester 
[5] , 4] , condensation of the vacancies and defects in solid 
Helium should result in a supersolid phase. It has how- 
ever, been argued l7| that the vacancies in '*He are too 
sparse to result in any meaningful phase transition. On 
the other hand, Leggett's proposal ,3] is tantamount to 
saying that supersolidity will result due to tunneling of 
real He atoms between neighbouring sites in the crys- 
tal. According to estimates based on his predictions, the 
supersolid fraction should be ~ 10~^ (Leggett's own es- 
timate was ~ 3 X 10^*) and the transition temperature 
< 1 mK. 

The results of the experiment |2] are, however, quite 
different. It is found that the transition occurs at 175 
mK and the supersolid fraction is 5 — 25 x 10"'^. We 
present here an argument that explains the emergence 
of the supersolid phase in He with the right parameters. 
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FIG. 1: Schematic Phase diagram of ^He. 



The theory, developed in _8] is peculiar to solid *He. It 
is a model to treat the local motion of atoms in solid 
^He. These modes are suugested to play a role in mass 
transport usually attributed to point defects (vacancies). 
The theory has also been successful in resolving the con- 
troversy regarding the contribution of point defects to 
the specific heat of solid ^He. From that we are also able 
to conclude that the helium in Vycor glass is essentially 
in a bcc phase and not in an hep phase. We also make 
contact with the concept of off-diagonal long range order 
(ODLRO) that is a necessary and sufficient condition for 
superfluidity in helium l lll |. 



We next review the essential points in the theoretical 
description of solid ^He as propounded in Ref . Q • The 
usual assumption concerning atomic motion in a solid is 
that the atom resides at a local minimum of the inter- 
atomic potential. In the quantum mechanical description 
however there is always a residual energy Ez even at zero 
temperature (zero point motion) which is well estimated 
by E, « -^, where m is the mass of the atom and a 
is the localization of the ground state. This makes itself 
most evident at the lowest densities. For Helium, Ez is 
so large that it remains a liquid under its vapour pressure 
down to absolute zero. Substantial pressure is needed to 
solidify it, as seen from Fig 1. At low temperatures, the 
minimal pressure necessary is 25 bar for ^He and some- 
what more so for ■^He. Much more pressure (40 bar) is 
needed to solidify helium in confined geometies [9|. It is 
perhaps the only solid whose phase diagram is so strongly 
influenced by quantum mechanics. The subsequent dis- 
course is valid only for low density solids (^--^ 21 cm^ molar 
volume). 

Within a small distance Ax around the minimum, the 
interatomic potential V{x) can be Taylor expanded as 



V{x) ^Vo + 



d^v 



dx^ 



(Ax) 



(1) 



where Vq is the minimum value and xq is the minima 
in the potential. However, for the solid helium atom 
potential there is a problem: the second derivative at the 
centre of the well is negative. A related problem in the 
conventional description is that the zero point motion is 
too large for the potential to be considered stationary. 

A key assumption of the harmonic treatment is that 
the nucleus is a point-like object as compared to the 
electron cloud. This approximation (Born-Oppenheimer 
approximation) is not valid when the potential is anhar- 
monic. In a quasi-classical picture, the nucleus will not 
remain at the centre due to the forces from the anhar- 
monic potential. This relative displacement gives rise 
to a dipole moment. Dipolar fluctuations of an atom 
in free space are random and isotropic. In a solid, the 
motion-induced dipole moments have preferred orienta- 
tions. Over time, the dipole-dipole interaction between 
these dipole moments would average to zero if the oscilla- 
tion frequency and phase of each of the atoms is random. 
If however, the zero-point motions are correlated, then 
the time averaged interaction is non-zero and can in cer- 
tain situations lower the energy of the solid. The lowest 
interaction energy is evidently when all the dipoles in 
a lattice are oscillating in phase. Since the direction of 
the dipole shows the instantaneous direction of the mo- 
tion, such a state is a translation of the system and hence 
unphysical. Looking for other symmetric arrangements, 
one finds two " antiferroelectric" configurations along the 
symmetry axes of the crystal with individual dipole mo- 
ments oriented along the (001) direction, as depicted in 




FIG. 2: The two "antiferroelectric" arrangements in the bee 
phase. Adapted from ||y|. 



Fig. 2 in the bcc lattice Il6| These have a total dipole 
moment of zero. Similarly, correlated motion can occur 
along the (010) and (100) axes too. It should however 
be noted that only one of these will exist in a region of 
space in the crystal to preserve this long-range ordering. 
The dipoles introduced above are affected by lattice 
vibrations, viz, phonons. For a modulation along a di- 
rection k, the dipolar interaction energy is 



x(fc) = -H 






3cos^[/x- (ro -Ti)] - 1 



\ro 



^2-Kik-{rQ —Vi) 
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where /i is the induced dipole moment and r^ is the in- 
stantaneous position of the i th atom. In the bee geom- 
etry, the only phonon modes that modulate the energy 
along the (001) direction are L(OOl), r(lOO) andTi(llO). 
It can be shown from periodicity properties that the cou- 
pling of the local modes to the lattice excitations is only 
through the ri(llO) mode Q- Thus the only elemen- 
tary excitations to the dipole array would be in the (110) 
direction, in the form of the ri(llO) phonon. 

For further analysis one introduces the Hamiltonian for 
the local interactions which are treated as bosons. The 
net Hamiltonian written in swcond quantized notation is 



H = J2[eo + X{k)](blbk + l)+J2xik)ibW-k + bkb-k) 

k k 

(3) 
where &J,,6fc are the bosonic creation/annihilation oper- 
ators. This Hamiltonian can be diagonalised using the 
Bogoliubov transformation (3k = u{k)bk + v{k)b_j^. The 
functions u{k),v{k), given in J8| are not essential to our 
analysis. The energy spectrum of the diagonalised Hamil- 
tonian is 



E(fc) = ^eo[eo + 2Xik)]. 



(4) 



Since the mode is an acoustic phonon, it must be gapless 
as fc ^ and hence X{k — 0) — ea. Using this 



e(fc) = eoV2sin[A; • {{re, — r^)] = eov2sinfco 



(5) 



in which the last inequality is particular to the T'i(llO) 
mode. 

As has been indicated in y, this mode can result 
in "phonon-assisted" tunnelling of the actual atoms be- 
tween the lattice sites. In that case, the atoms are in a 
'free-flowing' state that can undergo Bose Einstein con- 
densation and subsequently lead to supersolidity. This is 
essentially the mechanism suggested by Leggett [SJ. The 
difference is that we no more use the Hartree approxima- 
tion to calculate the dispersion relation of these modes 
as given in p]; we rather account for the anharmonicity 
through local modes to obtain the actual coherent diffu- 
sion process. 

One can calculate the supersolid transition tempera- 
ture from the elementary theory of non-interacting Bose 
gases. The average number of particles (bosons) is given 
as 



A^ = (n,) 



.g(e)c^£ 

j-ie/5£ - 1 



(6) 



where g{e) is the density of states. For a dispersion rela- 
tion given as 



e(fc) = ej, sinfco 



(7) 



putting e'„ = \/2eo 



He) 



V 



{sin-\/e'o)' 



a^27:^e'„ ^1 - {e/e',Y 



(8) 



with V being the specific volume (~ 21cc/mole). The 
lattice spacing in solid helium is a = 2.6 A. It immediately 
follows that 



N ^ 
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Since A^ is a monotonically increasing function of z, the 
condensation takes place at z = 1. Also the Landau 
critical velocity is given by the minimum value of ^^, 
where p = hk. Thus 



Vc = nun 



. (({p)^ ae^sinxo 



p 



J h 



Xq 



(10) 



Xq being a non-trivial root of tanx = x. Finally, a judi- 
cial estimate of the superfluid fraction is given by 



Ps 



' /fc2 /_„2\-l 



eo/E, = e'o{h'/ma') 



(11) 



A value of eg ^ 30^K or eq ~ 30/ v^ = 21.32^K repro- 
duces the results of the experiment J2| to excellent agree- 
ment. This value is also in harmonious agreement with 
the data from neutron scattering experiments, eo '^ 20/iK 
presented in |3| . There the magnitude of the coherent dif- 
fusion in solid helium is estimated. On the other hand, 
with the experimental value of eo — 20/iK, for a super- 
solid fraction of 10 x lO^'^, we obtain from © a tran- 
sition temperature of 215 mK and from (|10|l a critical 
velocity of 251/ims~^. The transition temperature com- 
pares well with the experimentally determined value of 
175 mK. The critical velocity is within the low tempera- 
ture limit of 300/ims~^. The introduction of the channel 
predicted in J3| accounts for the atomic self-diffusion in a 
way that is in commensuration with other facts like the 
specific heat of solid helium [lOj. Both these processes 
however call for the long range correlations characteristic 
of a quantum solid. 

The existence of off-diagonal long range order 
(ODLRO) is a necessary and sufficient condition for the 
existence of superfluidity 112] ■ The charecterization of su- 
persolidity used in |2| is the most literal one, taken from 
superfluidity: the incomplete dragging by slowly rotat- 
ing walls within which the material is confined. This is 
self-evident from their experimental setup. The proof in 
|12| about the necessity and sufficiency of ODLRO for 
superfiuidit y c an be trivially extended to the case of su- 
persolidity [ijl- Hence we can conclude that a system is 
supersolid or superfluid if and only if it exhibits ODLRO. 
Which of these will occur is of course determined by the 



thermodynamic coordinates of the system in the phase 
space, i.e, its temperature, pressure and density. 

From the Hamiltonian analysis previously carried out 
|8j , one finds the ground state wave function of the local 
modes to be 



■ 175mK 0. 
We would like to thank M. de Llano for helpful discus- 
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where U}^ , V}~ are the coefficients of the Bogoliubov trans- 
formation made in diagonalising the hamiltonian in Eqn. 
|3| This, as is well-known, is the famous Bogoliubov wave- 
function and known to exhibit ODLRO llj. The Bose- 
Einstein condensation is the simplest form of ODLRO 
and we believe that the emergence of the supersolid 
phase in sohd ''He is nothing but the BEG of the ri(llO) 
phonon mode. Therefore the existence of ODLRO in the 
local mode ground state wave function reaffirms the role 
of ODLRO in supersolid phase transition. 

In the present article we have showm that the emer- 
gence of a supersolid phase is a natural consequence of 
the quantum nature of solid ^He. Although such a phase 
was proposed a long time ago, the precise physical mech- 
anism for such transition was only recently presented. 
Experimental efforts in the mean time continued. We 
can now conclude that the supersolid transition in solid 
"'He is due to a Bose-Einstein condensation of the actual 
mass transport process amongst the helium atoms at dif- 
ferent sites in the bcc phase. The dispersion relation for 
this process is given by that of the ri(llO) phonon mode 
(or any of the two other orthogonal ones). For helium 
solidifying in the hep phase, there is no long range cor- 
relation (not to be confused with ODLRO) between the 
atomic dipoles and hence no phonon mode analogous to 
the one dealt with in this paper. Thus ^He solidifying in 
the hep phase will not exhibit a supersolid transition. For 
the experiment of solidifying helium in Vycor,as in Ref. 
[21 it might be tough to determine directly into which 
phase, bcc or hep the helium crystallizes. However the 
amount of helium going into the bcc/hcp phase can be de- 
termined from the value the supersolid fraction. As per 
our estimates, almost all of it goes into the bcc phase. 
We also show that the system is an ideal candidate for 
exhibiting a supersolid transition as it has ODLRO. 

As to ^He, as similar behaviour may be expected to 
exist due to the existence of an identical mode [l^ . The 
arguments leading to this mode are however valid only at 
temperatures much higher than the magnetic interactions 
(T ^ ImK). ^He is a fermion and the superfluid phases 
in the liquid state are observed at around 2 mK and 
presumably due to the formation of Cooper pairs. The 
exact temperature of transition can be hard to calculate 
due to the interactions in ^He [la|- Nonetheless, if such a 
supersolid phase is to exist it must be at temperatures at 
least 2 orders of magnitude below those for ''He. Indeed 
pure ^He shows no such transition at temperatures of 



sions. 
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